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Abstract 

Some systematic general features of y-scaling structure functions, which are 
essentially independent of detailed dynamics, are pointed out. Their physical 
interpretation in terms of general characteristics, such as a mean field descrip- 
tion and nucleon-nucleon correlations, is given and their relationship to the 
momentum distributions illustrated. A new relativistic scaling variable is pro- 
posed which incorporates the momentum dependence of the excitation energy 
of the (A — 1) system, with the resulting scaling function being closely related 
to the longitudinal momentum distributions and being free from removal- 
energy scaling violating effects. 
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Inclusive quasi-elastic electron scattering is potentially a powerful method for measuring 
the momentum distribution of nucleons inside a nucleus. Non-relativistically, this is most 
succinctly made manifest by expressing the data in terms of the scaling variable y which, 
over a large kinematic range, can be identified as the longitudinal momentum of the struck 
nucleon, k\\ [TJ. At sufficiently large momentum transfers, q, the structure function, W(u, q 2 ), 
which represents the deviation of the cross-section from scattering from free nucleons, scales 
to a function of the single variable y according to qW(u, q 2 ) pa f(y) where v is the electron 
energy loss and q = |q|. Thus, in the scaling limit, qW approaches a function that effectively 
traces out the longitudinal momentum distribution of the nucleons: 

oo 

f( y ) = J n(A;||,k ± )d 2 k ± = 2vr J n(k)kdk (1) 

Ivl 

Here, n(k) (with k = |k|) is the conventional nucleon momentum distribution function 
normalized such that 

/oo 
dyf(y) = 1 (2) 
-oo 

Knowledge of f(y) can therefore be used to obtain n(k) by inverting Eq. ([!]): 
The above picture relies on the simple assumptions that: 

i) nucleon binding does not play a role in the scattering process; in reality, qW(v, q 2 ) 
is determined by the spectral function, P(k,E), which depends both on the removal 
energy (E), as well as on the momentum of the nucleons, through the relation |2[ 
qW(y, q 2 ) = F(y) = f(y) - B(y) where 

B(y) = 2ir J dE J P 1 (k,E) (4) 

^min \y\ 

Pi being that part of P(k, E) generated by ground state correlations (thus, in a mean 
field description or, for the case of 2 H, Pi = 0) ||; 



ii) final state interactions (FSI) are disregarded; if they are taken into account a direct 
relation between the experimentally measured qW(u,q 2 ) and the asymptotic scaling 
function F(y) holds only approximately. 

Over the past several years there have been vigorous theoretical and experimental efforts 
to explore y-scaling over a wide range of nuclei Q, using a relativistic scaling variable 
resulting from energy conservation implicit in the instant form of relativistic dynamics (see 
Eq. |6| below) and representing the longitudinal momentum of those nucleons which have 
the minimal value of E (recently, it has been shown || that, for the deuteron, scaling in 
this variable y(= y r of ref. ||) also follows from a relativistic light-front approach). In Ref. 
the asymptotic scaling function F(y) has been obtained by an elaborate extrapolation 
procedure of existing data aimed at removing (or, at least, minimising) the effects of FSI. 
The longitudinal momentum distribution f(y) has thereby been obtained by adding to F(y) 
the binding correction B(y) evaluated theoretically. Such a procedure is based upon two 
basic assumptions: (i) the FSI can be represented as a power series in 1/q, and die out 
at large q 2 , a conclusion which has been reached by various authors || ;(ii) the theoretical 
binding correction has to be applied to obtain f(y) 0. Both assumptions, which in principle 
could be questioned, affect the final form of f(y), and therefore need to be investigated. We 
will discuss the effects of FSI in a separate paper 0; here we address only point (ii). To 
begin with, let us assume that f(y) obtained in is correct and let us analyze it in detail 
to see whether it contradicts or agrees with current theoretical predictions. A systematic 
analysis, to be presented elsewhere |7[], exhibits the following general features of f(y) for 
nuclei with A < 56: 

i) /(0) decreases monotically with A, from ~ lO^MeV^ 1 when A = 2 to ~ 
SxlO^MeV- 1 for 56 Fe; moreover, for y ~ 0, f(y) ~ C x [a 2 + y 2 )-\ with a (Ci) 
ranging from ~ 45MeV (18MeV) for A = 2, to ~ UOMeV (59MeV) for A = 56. 

ii) For 50MeV < \y\ < 200MeV, F(y) ~ e~ a2y2 with a ranging from ~ 56 x lO^MeV- 1 
for A = 2, to ~ 45xlO~ 2 Mey- 1 for A = 56. 
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iii) For \y\ > 400MeV, f(y) ~ C 2 e~ %l , with C 2 ranging from 2.5 x lO^Me^" 1 for A = 2, 
to 6x 10 _4 Mel/ -1 for A = 56, and, most intriguingly, b = 6 x 10 -3 MeV _1 , independent 
of A. 

The following form for /(y) yields an excellent representation of these general features 
for all nuclei: 

2 2 

fiv) = ^f^r + (5) 

The first term (= /o) dominates the small y-behavior, whereas the second term (= f\) 
dominates large y. The systematics of the first term are determined by the small and 
intermediate momentum behaviours of the single particle wave function. For \y\ < a this can 
be straightforwardly understood in terms of a zero range approximation and is, therefore, 
insensitive to details of the microscopic dynamics, or of a specific model. The small k 
behavior of the single particle wave function is controlled by its separation energy, (Q = 
M + Ma-i — Ma = E min ) and is given by (k 2 + a 2 )~ l so a = (2/iQ)^, /x being the reduced 
mass of the nucleon. 

Before discussing the intermediate range it is instructive to consider first the large y- 
behavior. Perhaps the most intriguing phenomenological characteristic of the data is that 
f(y) falls off exponentially at large y with a similar slope parameter for all nuclei, including 
the deuteron. Since (i) b is almost the same for all nuclei including A = 2, i.e., f(y), 
at large y, appears to be simply the rescaled scaling function of the deuteron; and (ii) 
b(m 1. 18/771) 1/chd(~ 4.35/m), we conclude that the term C2e~ f, ' s/ ' is related to the 
short range part of the deuteron wave function and reflects the universal nature of iViV 
correlations in nuclei. The remaining parameters, C% and a, can be related to /(0) and the 
normalization condition, Eq. (Q). Once this is done, there are no adjustable parameters for 
different nuclei. The intermediate range is clearly sensitive to a, with the gaussian form 
being dictated by the shell model harmonic oscillator potential, modulated, however, by the 
correct \y\ < a behaviour, namely (y 2 + a 2 ) -1 , thereby ensuring the correct asymptotic wave 
function. As an example, Fig. 1 shows the longitudinal scaling functions f(y) for 2 H and 



A He extracted from the experimental data ( @), compared to Eq. The errors of f(y) for 
2 H are very small, whereas, at large values of \y\ they are appreciable for 4 He and heavier 
nuclei; they are due both to the lack of reliable experimental data at large values of q 2 , and 
to the necessity to correct the asymptotic scaling function for binding effects, and produce 
an error on C2 and b in Eq. ([5]), ranging from ~6%to~10%. A systematic analysis for a 
large body of nuclei exhibiting the same features as those shown in Fig. 1 will be presented 
elsewhere @. 

With these observations it is now possible to understand the normalization and evolution 
of f(y) with A. First note that Eq. (|) implies /(0) = \ f d 3 k^ = (l/2k) and so is mainly 
sensitive to small momenta. Now, typical mean momenta vary from around 50 MeV for 
the deuteron up to almost 300 MeV for nuclear matter. We can, therefore, immediately 
see why /(0) varies from around 10 for the deuteron to around 2-3 for heavy nuclei. More 
specifically, since C2 <C C\/a 2 and f\ falls off so rapidly with y, the normalization integral, 
Eq. (|2]), is dominated by small y, i.e., by f . This leads to /(0) « (^^a)^ 1 = (2n^Q)~ 1 ^ 2 
which gives an excellent fit to the A-dependence of /(0). Since f(y) is constrained by the 
sum rule, Eq. (0), whose normalization is independent of the nucleus, a decrease in /(0) 
as one changes the nucleus must be compensated for by a spreading of the curve for larger 
values of y. Thus, an understanding of f(y) for small y coupled with an approximately 
universal fall- off for large y, together with the constraint of the sum rule, leads to an almost 
model-independent understanding of the gross features of the data for all nuclei. 

To sum up, the "experimental" longitudinal momentum distribution can be thought 
of as the incoherent sum of a mean field shell-model contribution, (/ ), with the correct 
model-independent small 2/-behaviour built in, and a "universal" deuteron-like correlation 
contribution (fx). Thus, the momentum distribution, n(k), which is obtained from (|3]), is 
also a sum of two contributions: n — no + Tii. This allows a comparison with results from 
many body calculations in which uq and U\ have been separately calculated. Of particular 
relevance are not only the shapes of no and rix, but also their normalizations, Sam = 
/ n (i)(i 3 fc = J fo(i)dy which, theoretically, turn out to be, for 4 He, So ~ 0.8 and S\ ~ 0.2 || 
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whereas Eq. (|5]) yields So = 0.76 and Si = 0.24. A comparison between the momentum 
distributions obtained from y-scaling and the theoretical ones is shown in Fig. 2. As can be 
seen the n(k) compare very well with theoretical calculations. Our analysis confirms very 
well known properties of low energy nuclear physics (binding, etc), as well as some features 
of the momentum distributions predicted recently by various theoretical calculations |§ , @ , 
III , but still waiting for a firm experimental confirmation. To place our results on a more 



solid basis, it would be necessary however to reduce the errors on f(y) at large values of 
\y\, which necessitates experimental data at larger values of q 2 , as well as a reduction of 
the uncertainties related to the binding correction. Experimental data at higher Q 2 became 



recently available and their inclusion in the analysis reduces indeed the errors on f(y) 
0. Here we address the problem of making the extraction of f(y) as much independent as 
possible from theoretical binding corrections. This is accomplished by introducing another 
scaling variable. A generic ambiguity in defining a scaling variable is that there is no unique 
prescription for y, so that it is legitimate, in principle, to incorporate in the definition of 
y some physical dynamical effects by introducing proper effective parameters. The usual 
scaling variable y is obtained from relativistic energy conservation 

v + Ma = [(Ma-i + E* A _,) 2 + k 2 ] 1/2 + [M 2 + (k + q) 2 ] 1/2 (6) 

by setting k — y, ^ = 1, and, most importantly, the excitation energy, E\_ x = 0; thus, y 
represents the nucleon longitudinal momentum of a nucleon having the minimum value of the 
removal energy (E = E min , E* A _ X = 0). The minimum value of the nucleon momentum when 
q -> oo, becomes k min (y,E) = \y- (E-E min )\. Only when E = E min does k min (y,E) = \y\, 
in which case the binding correction B = and F(y) = f(y). However, the final spectator 
(A — 1) system can be left in all possible excited states, including the continuum, so, in 
general, E* A _ X ^ and E > E min , so B(y) ^ 0, and F(y) ^ f{y)- Thus, it is the dependence 
of k m in on E* A _ X that gives rise to the binding effect, i.e. to the relation F(y) ^ f(y)- This is 
an unavoidable defect of the usual approach to scaling; as a matter of fact, the longitudinal 
momentum is very different for weakly bound, shell model nucleons (for which E A _ 1 ~ 



— 20MeV) and strongly bound, correlated nucleons (for which E* A _ X ~ 50 — 200MeV), so 
that at large values of \y\ the scaling function is not related to the longitunal momentum 
of those nucleons (the strongly bound, correlated ones) whose contributions almost entirely 
exhaust the behaviour of the scaling function. In order to establish a global link between 
experimental data and longitudinal momentum components, one has to conceive a scaling 
variable which could equally well represent longitudinal momenta of both weakly bound and 
strongly bound nucleons. An attempt in such a direction has been made in the past by |14 



and, recently, by us |H| with some minor differences with respect to Ref. attempt. This 



was based upon taking literally the two-nucleon correlation model according to which the 
large k and E behaviours of the Spectral Function are governed by configurations in which 
the high momentum of a correlated nucleon (1, say) is almost entirely balanced by another 
nucleon (2, say). Within such a picture, one obiously has E A-l = -jzjqJi^ 2 and the average 
excitation energy for a given value of k is < E^^ik) >= jrE^IM- By replacing E* A-l in 
Eq. (||) with < _E^_ 1 (fc) >, the deuteron-like scaling variable y 2 introduced in our previous 
paper |I3|(see also [TJ]], where a similar scaling variable was first introduced) is obtained, 



representing the scaling variable pertaining to a "deuteron" with mass M = 2M — EH', 

(2) 

where E\ h = Ma-2 + 2M — Ma- Such a scaling variable, however, has the unpleasant 
feature that the effect of the deuteron- like correlations are overestimated at low values of 
U2 and , as a result, the correct shell-model picture provided by the usual variable y is lost. 



In a more refined model, the CM motion of the pair is taken into account: one has [|TT 

which shows that the excitation energy of the residual nucleus depends both upon k and 
K CM ; although using Eq. (|7|) the situation is improved, the description ot low values of y is 
not a satisfactory one. In this paper we assume that the nucleus is described by a realistic 



spectral function as provided by few- and many-body calculations |12|], [|13j and implement 
such a realistic description into the definition of the scaling function. If the expectation 
value of Eq. (0) is evaluated with realistic spectral functions, for the three-body system [12| 



and nuclear matter ]T3[ and the model spectral function of JTT| one obtaines ( II7IJ11 

1 



CA 



2M 



(8) 



Here, b A and ca, resulting from the CM motion of the pair, have values ranging from 
17MeV to 43MeV and 3.41 x 10 _1 to 1.66 x 10 _1 , for 3 He and Nuclear Matter, respectively. 
Placing Eq.(|8|) in Eq.(^|) and subtracting the value of the average removal energy < E > 
to counterbalance the value (§) at low values of y, a new scaling variable is obtained. This 
effectively takes into account the fc-dependence of the excitation energy of the residual 
(A — 1) system, both at low and high values of y. This is in contrast to the usual scaling 
variable, which completely disregards E* A _ Xl and the scaling variable y2, which overestimate 
the effects of deuteron-like correlations at low values of y. In the kinematical region of 
existing experimental data this new , global scaling variable, yew-, that we have obtained, 
reads as follows 

1/2 

^l/A IV1 — VVA 
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yew 



4W A ' 



Here, v A = v + M, M = (2A - 3)M/(A - 1) - E$ - {b A + 2M 2 c A - < E >), q = q — c A v A 
and Wl = z/ A 2 — q 2 = M 2 + 2z/M - Q 2 . For the deuteron E\_ x = 0, so y cw -»• V = I - ?/2 + 
[g 2 /4 - (WM 2 - W d 4 )/W d 2 ] 1/2 | with u d = u + M d and W 2 , = u d 2 - f = M d 2 + 2uM d - Q 2 . 
For small values of ycw-i sucn th & t (4zf ^^a — ca-^m) ^< E >, the usual variable, 



(9) 



representing the longitudinal momentum of a weakly bound nucleon is recovered ||21|| . Thus 
yew interpolates between the correlation and the single particle regions. More importantly, 
however, since k min (q,u,E) ~ \y C w\, B(y C w) — 0, F(ycw) - f{vcw)- One would therefore 
expect from our above analysis, the same behaviour of f(ycw) & t high values of yew f° r both 
the deuteron and complex nuclei (unlike what happens with the usual scaling function F(y)), 
and the same shell-model behaviour at low values of y, as predicted by the usual scaling 
variable. This is, indeed, the CclSC, clS exhibited in Figs. 3 and 4, where the direct, global, 
and independent of A link between the scaling function F(q, yew) and the longitudinal 
momentum distributions is manifest . 



S 



We can summarise our conclusions as follows: 

i) The general universal features of the y-scaling function have been identified and in- 
terpreted in terms of a model-independent zero-range contribution, a "universal" 2- 
nucleon correlation contribution and a mean field (shell-model) contribution. The 
shape and evolution of the curve have been understood both quantitatively and qual- 
itatively on general grounds. 

ii) A global relativistic scaling variable which, unlike all previously proposed variables, 
incorporates the average mean field excitation energy of the (A — 1) system, as well as 
the excitation energy produced by deuteron-like correlated pairs and by their center- 
of-mass motion in the nucleus, has been defined. Such a variable, thanks to these new 
features, allows one to establish a more direct link between the scaling function and 
the longitudinal momentum distributions at any value of y. Thus, using this variable, 
it would be possible in principle to obtain the longitudinal momentum distributions di- 
rectly without introducing theoretical binding corrections. Of course the usual variable 
has the advantage of being defined in terms of a well-defined experimental quantity, 
the minimum value of the removal energy E min , whereas yew incorporates the exci- 
tation energy of (A — 1) by a theoretical prediction. However, given the fact that in 
7* — Nucleus scattering the virtual photon couples to nucleons having different values 
of the removal energy, suggests that the removal energy has to be taken into account 
in the definition of the scaling variable. Our Fig. 3 shows indeed that removal energy 
effects are very important. Since these are a source of scaling violation, the other 
source being the FSI, it seems reasonable to incorporate the binding effects into the 
definition of y so as to ascribe the remaining scaling violation to FSI. Moreover the 
plot of the data in terms of yew has a very clear cut meaning: F(ycw, <? 2 ) represents 
the scaling function at a given q 2 and for a longitudinal momentum of a nucleon having 
removal energy E { y cw ) = E$ + j-^^ycw 2 + bA-c A ^ycw As already pointed out, 
there is no unique prescription for defining a y-scaling variable, and various variables 
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are currently being used [VJ§. Ultimately, a valid criterion for a scaling variable, is to 



produce scaling; the usual y has been shown to produce scaling but, at large values of 
\y\, in the ultra asymptotic limit ||, ||18|| , that is a limit which is not reached by 



present experimental data. In Ref. 0, it is shown that the new scaling variable not 
only produces precocious scaling to the longitudinal momentum distribution even at 



the largest value of \y\ recently reached by the new TJLAB experimental data ||19 
but also that the FSI on the scaling function F(ycw, <? 2 ) are very similar to the ones 
acting in the deuteron. Therefore, in terms of this variable the data seem to support 
the idea that the large yew behaviour in all nuclei is essentially nothing but a rescaled 
version of the deuteron (including, perhaps, also some effects from deuteron-like FSI, 
due to the constant value of the nucleon-nucleon cross section in the region 2 < Q 2 < 6 
GeV 2 , as stressed in p0| ). Such a conclusion cannot be reached by analysing the data 
in terms of the old scaling variable, which mixes up scaling violation due to removal 
energy effects and FSI. 

We would like to thank Avraham Rinat and Dino Faralli for useful comments and 
discussions. 
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FIGURES 

FIG. 1. The "experimental" longitudinal momentum distributions of 2 H and 4 He ob- 
tained using the results of (2| compared with Eq. (|5|) with a = 45MeV, C\ = 18MeV, 
a = 56 x lQ-iMeV' 1 , C 2 = 2.5 x lO^MeV' 1 , and b = 6 x lO^MeV -1 , for 2 #, and 
a = 167MeV, Ci = 106MeV, a = 68.5 x W^MeV' 1 , C 2 = (6 ± 0.6) x lO^MeF" 1 , and 
6 = (6 ± 0.6) x Kr 3 MeV _1 , for 4 #e 

FIG. 2. The nucleon momentum distribution for 2 H from Eq.(||), compared with the one 
obtained from the AV14 interaction; the same as Fig. 2a but for 4 He 

FIG. 3. The experimental scaling function F(q, y) for 2 H, 4 H and 56 Fe compared with the lon- 
gitudinal momentum distributions f(y) given by Eq.Q (dot-dash- 2 .?/; short-dash- 4 i/; full- 56 .Fe). 
The scaling variable is the usual one, i.e. the one obtained from Eq.(|5|) placing E^_ 1 = 

FIG. 4. The same as in Fig.3, but with the variable , yew, obtained using Eq.(||) 
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Fig. |T|. Ciofi-West Y-scaling 
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Fig. £|. Ciofi-West Y-scaling 
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